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ABSTRACT 

We present a study of the stcUar populations and morphologies of galaxies in the core 
of the galaxy cluster Abell 2218. Integral field spectroscopy (IFS) observations were 
performed using PMAS in the PPAK mode covering a field-of-view of ~74" x64" cen- 
tred on the core of the cluster, in order to obtain spectroscopy of an unbiased flux 
limited sample of cluster galaxies. Forty-three objects were detected in the IFS data, 
31 of them with enough signal-to-noise to derive the redshift, all of them brighter 
than I <21.5 mag. Twenty-eight are at the redshift of the cluster (17 with previously 
unknown redshift). Individual spectra of the cluster members were extracted and com- 
pared with single stellar population models to derive the luminosity-weighted ages and 
metallicities. In addition, deep HST/ACS F475W, F555W, F625W and F850LP-band 
images centred on the cluster core were obtained from the HST archive (zum ~28 
mag). A detailed morphological analysis of all the galaxies within the field-of-view of 
these images down to znm <22.5 mag was performed classifying them as late-type, 
intermediate and early-type, on the basis of their Sersic indices. The literature was 
scanned to look for spectroscopically confirmed cluster members located within the 
field-of-view of the ACS image. The final sample of 59 galaxies comprises our reported 
sample of 28 galaxies in the core, and 31 additional galaxies in the outer regions. In 
addition, multiband broad-band photometry was extracted from the literature for all 
objects. 

The morphologically segregated color-magnitude diagram shows that the early- 
type galaxies cover the range of brighter and redder colors (the so-called "red se- 
quence"). A large fraction of spiral galaxies (~50%) is found, most of them fainter 
than the limit of previous studies. They cover a wide range in colors, from blue colors 
typical of Butcher-Oemler galaxies to red colors similar to those of early-type galaxies. 
This result indicates that early-type galaxies are more massive and have older stellar 
populations, while late-type galaxies are less massive and have a wider range of stellar 
populations. The distribution of luminosity weighted ages as a function of metallicities 
and luminosity weighted masses, and the distribution of Sersic indices as a function 
of the luminosity weighted masses confirm these results. They in fact agree with a 
proposed two-step scenario for the evolution of galaxies in clusters, where the star 
formation is quenched first in the infalling spirals, after which a morphological trans- 
formation follows that requires larger time-scales. This scenario naturally explains the 
population faint late-type galaxies with old stellar populations observed in this cluster. 
In addition, an extremely blue merging galaxy system is found at the core, with the 
nominal redshift of the cluster. 

Key words: galaxies: clusters: individual: Abell 2218 - galaxies: elliptical and lentic- 
ular, cD - galaxies: spiral - galaxies: irregular - galaxies: fundamental parameters - 
galaxies: stellar content - techniques: spectroscopic 
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1 INTRODUCTION 

Galaxy clusters have been used for decades to study the 
evolution of galaxies. Being tracers of the largest density 
enhancements in the Universe, they are considered the re- 
gions where galaxies formed first. At low redshift they are 
dominated by large early- type galaxies (E/SO), which col- 
ors are consistent with a bulk formation of their stars at 
high redshift. Classically, it was thought that the evolution 
of this dominant population through different cosmologi- 
cal epochs was largely passive without signatures of im- 
portant deviations (ie., periods of more recent star for- 
mation). Howeve r, this sim p le pic ture has change d in the 
last d ecades. Butcher et al.l (|l983l ) and [Butch er fc Oemlerl 
l| 19841 1 showed an increase in the fraction of blue galaxies 
in clusters with redshift. High-resolution ground-based and 
HST images have shown that they correspond to galaxies 
with late-type morphologies, often distur bed or with close 
companions (e.g. Tho mpson 1986 , 1988 : .Couch et al.lll994 
[Dressier et al.l Il994l . Il997l ) . One observes a morphological 
evolution of the galaxies in clusters, in the sense that the 
fraction of E-type galaxies remains constant as a function 
of different cosmological epochs, while the fraction of SO's 
incre ases at the expens e of L-types towards lower redshift 
(e.g. iFasano et ahlbOOOl ). 

The presence in clusters of a significant population of 
post-starburst galaxies known as "E-l-A" or "k-t- a" galaxies 
IIDressler fc Gunnlll983l : ICouch fc Sharpies! 1 198'iD . indicates 
that star formation was active in the recent past a nd sharply 
halted at s ome point during the last 1-1.5 Gyr l|Poggiantil 
12001 [ 20041 1 . Most of these galaxies are late type, and they 
show a radial distribution within the clusters that is inter- 
mediate bet ween the passive and the emission-line galaxy 
populations (jPoggianti et al.lll999l) . These results lead to 
the hypothesis of a two-step evolution of galaxies. A spiral 
galaxy is captured by the cluster and, during its infall, suf- 
fers from gas-st ripping due to the tidal forces involved (e.g. 
IVerheiie"nlll997h . Provoked by the lack of gas, the star for- 
mation is halted, and the galaxy evolves passively, being ob- 
served as a "k-|-a" post-starburst. An initial short starburst 

J irocess, previou s to the gas-stripping, may be also present 
Poggiantill2004l '). After this the galaxy evolves morphologi- 
cally over larger times scales. The disk dims due to the lack 
of a new generation of stars, that otherwise will dominate 
its luminosity, and the galaxy transforms itself into a lentic- 
ular galaxy. Dry mergers between lenticula r and spheroida l 
galaxies build up new large spheroidals (e.g. lBell et al.ll2006l ) 
The proposed scenario suggests that a fraction of 
the observed SO's were previously star-forming spirals. For 
these galaxies the star formation was truncated at some 
time between 2 and 5 Gyr ago and it is expected that 
thei r luminosity faded by ~0.5-1.5 mag and become red- 
der l|Poggianti et al.ll200lf ). populating the red sequence at 
magnitudes fain ter than the brighest 1 or 2 magnitudes 
l|Poggiantil l200^). The same process may affect dwarf galax- 
ies captured by the cluster, which will suffer more strongly 
the effects of tidal interaction and galaxy harr asment, which 
induces a short b ut intense starburst ( Moore et al.l 1 19961 : 
iRakos et al.ll200ll ). Most of the field dwarf galaxies are gas- 
rich rich , blue galaxies with late-type morpholo gies at any 
redshift l|Cair6s et al.ll200ll :lB arazza et al.ll2006D . Once cap- 
tured by the cluster, consumed or stripped of its gas, it is 



Table 1. Summary of the HST/ACS images 
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expected that their luminosity fades and they become red- 
der too, populating the faint-end of the red sequence. In 
both cases ones expects to find a population of red late- 
type galaxies, fainter than the original early-type galax- 
ies in the cluster. However, most of the previous studies 
have found tha t late-type galaxies are mostly blue (e.g. 
lThompsonlll986h . 

To test this hypothesis, complete spectroscopic surveys 
of cluster galaxies and high-resolution deep images for a 
proper morphological analysis of the faintest possible tar- 
gets are required. Complete spectroscopic surveys are dif- 
ficult and time consuming, due to the requirement of de- 
blending and placing objects in the spectrograph multislits 
or masks. In particular, the blending of different components 
in the core of the clusters is a problem, impossible to han- 
dle with classical slit spectroscopy. Due to that, most of the 
spectroscopic surveys are based on pre-selection of the clus- 
ter candidates, most of them based on color or morpholog- 
ical selections (e.g . [Dressier et al. 199 9; Balogh o t al. 1999|; 
iKelson et al]|200ol : Iziegler et al.ll200ll : iFritz et"aLir2'005 l. Al- 
ternatively, they are focused on galaxies in the less crowded 
outer regions. These selections may affect substantially the 
conclusions. For these reasons we have started an obser- 
vational program aimed at obtaining Integral Field Spec- 
troscopy (IFS) of the core of some of the richest galaxy clus- 
ters in the local universe, including Abell 1689 (ESO 73. A- 
0605) and AbeU 2218. For both of them deep HST/ACS 
images are available. A simi lar approach has been used be- 
fore bv lCovone et al.l (|2006l ). to study the dynamical mass 
of low redshift galaxy clusters. 

Ab ell 2218 is one of t he richest clusters in the Abell cat- 
alogue ijAbell et al.lll989l ). with a richness class of 4. It has 
a re dshift of z ^O-l/' ^'^ '^ a velocity dispersion of 1370 km 
s"^ (iKristian et al.lll978l : He Borgne et al.lll992h . The clus- 
ter sh ows a set of arc-lenses. Detailed high-resolution X-ray 
maps (|McHardv et al.lll990l ) and mass-concentration stud- 
ies based on the properties of the gravitational lenses have 
shown that the cluster has two density peaks, the strongest 
of them dominated by a large cD galaxy. Several observa- 
tional programs have provided us with a large dataset of well 
deblended slit spectra for the galaxies in the outer part of the 
cluster (jZiegler et al.., 20011) . as well as very good mul t iband 
ground-based arid HST imaging (iLe Borgne et al ] Il992l : 
ISmail et al.ll200ll : IZiegler et aLll200lh . We focused our sur- 
vey on the central arcmin region of the cluster, around the 
cD galaxy. 

The outline of this article is as follows. Section 2 de- 
scribes the observational dataset, including a description of 
the IFS data and the HST/ACS images. Section 3 includes 
the analysis of the dataset and presents the results. Here, 
Sections 3.1 and 3.2 describe the automatic morphological 
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Figure 1. Average night-sky emission spectrum derived by com- 
bining the spectra obtained through the 36 sky-fibers for the dif- 
ferent pointings of each night. 



analysis that we performed as opposed to classical methods. 
Section 3.3 shows the method used to deblend and extract 
the integrated spectra of the different detected objects from 
the IFS data, with the redshift determination explained in 
Section 3.4. The analysis of the morphologically segregated 
color-magnitude diagram is included in Section 3.5. The lu- 
minosity weighted ages and metallicities of the galaxies, and 
their relation to their masses and morphologies are discussed 
in Section 3.6 and 3.7. The results are discussed in Section 
4, and the conclusions are listed in Section 5. Throughout 
this article we have used a A-cosmology, with Hq—IQ km s^^ 
Mpc"\ J7m=0.3 and Q,a^Q.7. 



2 OBSERVATIONAL DATASET 

2.1 IFS observations and data reduction 

Observations were carried out on June the 30th and July 
the 6th, 2005, at the 3.5m telescope at Calar Alto Observa- 
tory with the Pots dam MultiAperture Spectrograph, PMAS , 
dRoth et al.ll2005l 'l in the PPAK mode (|Verheiien et al.ll2004 
iKelz et alT 20061 ). The V300 grating was used, covering a 
wavelength range between 4650-8000 A with a resolution of 
~10 A FWHM and a sampling of 3.2 A/pixel. The PPAK 
fiber bundle consists of 382 fibers of 2.7" diameter each one. 
Of them, 331 fibers (the science fibers) are concentrated in a 
single hexagonal bundle covering a field-of-view of 72" x 64" , 
with a filling factor of ~65%. The sky is sampled by 36 ad- 
ditional fibers, distributed in 6 bundles of 6 fibers each one, 
located following a circular distribution at ~90" from the 
center and at the edges of the central hexagon. The sky- 
fibers are distributed between the science ones in a pseudo- 
slit, in order to have a good characterization of the sky. The 
remaining 15 fibers are used for calibration purposes, as de- 
scribed below. 

For each night, three dithered exposures of 1 hour 
each were taken, following a pattern of ARA= -f 1.56" and 
ADEC= ±0.78", in order to have a complete coverage of the 
science field-of-view and to increase the spatial resolution. 
Each pointing was divided into shorter individual frames for 
a proper cosmic-ray removal. A total of 6 hours integration 



time was achieved after combining the exposures of the two 
nights (3 hours per night). 

Data reducti on w as perf ormed using R3D 

l|Sanchez fc Cardiell l2005l : ISanched [20 06) , in combina- 
tion with the IRAF package (|Todvl Il993l fl and E3D 
(jSanchez et al.l l2004al '). The reduction consists of the 
standard steps for fiber-based integral-field spectroscopy. 
A master bias frame was created by averaging all the bias 
frames observed during the night and subtracted from the 
science frames. The location of the spectra in the CCD 
was determined using a continuum illuminated exposure 
taken before the science exposures. Each spectrum was 
extracted from the science frames by coadding the fiux in 
5 pixels in the perpendicular direction of e ach fiber, and 
stored in a row-staked-spectrum file RSS (jSanchez et al.l 
l2004al ). Wavelength calibration was performed using a He 
lamp exposure taken at the beginning of the night, and 
corrected for distortions using ThAr exposures obtained 
simultaneously with the science exposures through the 
calibration fibers (indicated above). The final accuracy of 
the wavelength calibration was estimated by comparing 
the central wavelength of the sky-emission lines with their 
nominal values, finding an agreement within a range of ~0.3 
A. Differences in the fiber-to-fiber transmission throughput 
were corrected for by comparing the wavelength-calibrated 
RSS science frames with the corresponding continuum 
illuminated images. Flux calibration was performed using 
standard calibration stars observed throught the night 
(BD33-I-2642 and Hz44). The uncalibrated spectra of the 
stars extracted from the reduced frames were then compared 
with publically available calibrated spectra, determining a 
fiux calibration correction that was applied to the science 
frames, once corrected by the difi'erences in exposure time. 
The accuracy of the spectrophotometric calibration was 
~0.15 mag across the wavelength range covered by our 
dataset, as determined by comparing with high accurate 
broad-band photometry , as we will discuss below. A night 
sky-emission spectrum was obtained for each pointing by 
combining the 36 spectra obtained through the sky-fibers, 
and subtracted from the 331 science-fibers spectra. Figure 
[1] shows the average night-sky spectrum obtained after 
combining all of them. Several strong night-sky emission 
lines contaminate our spectra, in particular in the reddest 
part (A >7200A). 

The three dithered exposures were then combined, cre- 
ating a single RSS frame with 993 science spectra and 
the corresponding combination of the position table. A 
final datacube with 1" /pixel sampling was created for 
each nig ht by interpolating th e combined RSS frame us- 
ing E3D (|Sanchez et al.l [2"004al '). This corresponds to 2.85 
kpc/pixel at the redshift of our target. Finally, the two 
datacubes obtained for each night were recentered at each 
wavelength, correcting for differential atmospheric refrac- 
tion, and combined using IRAF tasks. The final datacube 
comprises 4828 individual spectra, covering a field-of-view 
of ~71"x68". 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 2. Grayscale image of the combined F850LP-band image, together with labels indicating the location of the galaxies studied 
by Le Borgne et al. (1992) (orange stars), Ziegler et al. (2001) (blue squares) and Small et al. (2001) (red crosses). The location of the 
confirmed cluster members sampled by our IFS data are also indicated (green circles). North is up and East left. 



2.2 HST/ACS images and literature data 

We obtained publically available HST/ACS images of the 
cluster in the F475W, F555W, F625W, F775W and F850LP- 
bands from the Hubble Space Telescope archive. Table [T] 
summarizes their main characteristics, including the cen- 
tral wa velength and width of their bands (obtained from 
ISirianni et al.ll2005l ). They basically correspond to images in 
the the classical B, V , R, I and z bands. The images were 
calibrated to the Vega system using the corre sponding mag- 
nitude zeropoints and prescriptions listed in ISirianni et al.l 
l|2005l ). Table [T] also lists the number of individual frames ac- 
quired, each one with a different orientation and a slightly 
different pointing, due to the observing strategy of the 
HST/ACS camera (in order to cover the gap in-between the 
two different CCDs). The acquired frames were already pro- 



cessed by the HST/ACS reduction pipeline, which includes 
bias and dark current correction, flat-fielding, quantum effi- 
ciency correction, field-distorsion correction, the derivation 
of the astrometric-solution and flux calibration. We correct 
each individual from for cosmic rays, realigned to the stan- 
dard North-East and finally combined into a single image 
per filter after a sky-background determination and sub- 
traction. The final exposure time and limiting magnitude 
of each image, estimated by looking for the faintest object 
detected at ~15cr over the background, are listed in Table 
[T] All the images are much deeper than our IFS data, as we 
will see below in detail. In particular this is the case for the 
F850LP-band image. This image is much deeper than the 
rest, being as de ep as the corresponding image of UDF (e.g. 
ICoe et aLlBooel ), and therefore one of the deepest images 
ever taken of this cluster. The cluster core is centred on one 
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Figure 3. left-panel: Three-color image created by coadding the flux of the flnal datacube through three broad-bands corresponding 
approximately to B, R and /. right-panel: Similar image created using HST/ACS data described in the text. An electronical version of 
the full size image can be found in the webpa.ge \http://www.caha.es/sanchez/abell2218/BRI^big.jpg\ 



of the ACS chips for this particular image, and the field- 
of-view of our IFS data is completely covered by it. This 
high resolution ACS F850LP-band image will be used for a 
detailed morphological analysis. 

We also collected broad-band photometry and redshift 
infor mation of the cluster g alaxies available in the litera- 
ture. iLe Borgne et al.1 (|l992l ) performed a photometric and 
spectroscopic survey of the galaxies in a field of 4'x4' 
centred on Abell 2218. This survey comprises broad-band 
photometry in 5 bands {B,g,r,i and z filters) for 729 ob- 
jects, and redshift determination for 66 of them, 51 clus- 
ter members. They calibrated their photometric data us- 
ing two differe nt systems, the Johnson one for the _B-band 
l| Johnson fc M orga n 1951^). and the T huan-Gunn for the 
g,r,i and z-band (|Thuan fc Gunnll 19761 ). We did not perfom 
any recalibration of their data to a common photometric sys- 
tem. These data were initially used to compare with the cor- 
responding intensity of the individual spectra of the objects 
detected in our IFU data, and for doing so we just needed 
to use th e correspoding Johns on and Thuan-Gunn effective 
flux (e.g iFukugita et al.|[l995l ). The accuracy of their pho- 
tometry ranges from 0.1 mag to 0.5 mag from the brighter 
{B ~20 mag) to the fainter {B --^ 26 mag) objects. The 
spectroscopic data lacks the required quality needed to per- 
form an accurate analysis of the age-metallicity indicators, 
like the equivalent widths of H/3 and Mgb, but they provide 
high quality redshift determinations. 

IZiegler et all (|200lh reported on multislit moderate- 
resolution spectroscopic observations of 48 color-selected 
early-type members of the cluster. They also present high- 
quality multiband photometry in 4 bands {U, B, V and /), 
derived from groud-based observations, and a morphologi- 
cal analysis of 19 of the 48 galaxies based on H ST/WFPC2 
observations (F702W-band, ~ ii-band). Finallv. lSmail et aP 
l|200ll ) present a photometric study of the ages and metallic- 
ities of 81 early-type cluster candidates (33 with confirmed 



redshifts), based on ground-based Jf-band and HST/WFPC 
multiband images, in the F450W, F606W and F814W- 
bands. In both cases, their photometric data were calibrated 
in the Vega-based system. Many other studie s have targeted 
this cluster since its redshift confirmation by Kristi an et al. 
(1978), and the early photometric studies bv iButcher et al. 
(jl983r ). However, in this paper we only use the photometric 
and spectroscopic data of the articles listed above. 

Figure [5] shows a greyscale of the combined HST/ACS 
F850LP-band image, together with labels indicating the lo- 
cation of the galaxies studied by Le Borgne et al. (1992) (or- 
ange stars), Ziegler et al. (2001) (blue squares) and Small 
et al. (2001) (red crosses). The location of the confirmed 
cluster members sampled by our IFS data are also indicated 
(green circles). As we will see below, most of them have not 
been previously targeted by any spectroscopic survey. 

2.3 The depth of our IFS data 

To obtain a rough estimate of the depth of our IFS data, we 
compared them with available deep optical images. Figure[3] 
shows a three-color image created by coadding the flux of the 
final datacube through three broad-bands corresponding ap- 
proximately to B, R and I (left panel), together with a sim- 
ilar image created using the previously descrived HST/ACS 
dataset. It is interesting to note the similarities between the 
two images, despite the fact that our datacube has a final 
FWHM for point-like sources of ~2.5" , much larger than the 
one from the ACS. Based on the /-band magnitude of the 
faintest detected object, we roughly estimate the detection 
limit for our IFS data Inm ~22.4 mag. However, we only 
get high enough quality spectra for the proposed study for 
objects brighter than I<21.5 mag (signal-to-noise ratio >5a 
per spectral pixel), as we will see below. This limit is far be- 
yond previous s pectroscopic studies on this cluster, around 
him ~18.7 mag l|Ziegler et al.ll2001^ . 
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Table 2. Summary of the photometric analysis. 
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3 ANALYSIS AND RESULTS 
3.1 Photometric analysis 

We used the previously described HST/ACS images to ob- 
tain homogeneus multiband photometry for aU the detected 



object in field of Abell 2218. This homogeneus photomet- 
ric dataset is fundamental to estimate the accuracy of the 
the spectrophotometric calibration of our IFS data and to 
study the stellar pop ulations of the gala xies in the cluster. 
We use SExtractor (|Bertin et al.l Il996l) to detect objects 
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Figure 4. Results from the simulations. Solid-dots show the re- 
covered Scrsic indices of the simulated galaxies as a function of 
the magnitude difference between the bulge and total components 
(Amag). Overplotted countours show the density distribution of 
the plotted points. 

in the images up to 15cr over the sky (20 connected pixels 
with a signal-to-noise of at least 2a each one). This limit 
is very relaxed, and most probably we are losing many real 
faint objects by using it. However, we are not interested 
in the faintest targets, most probably high-z objects or ar- 
clenses (see Fig 1), but only in the cluster members. For each 
detected target we obtained aperture photometry in each 
band, using the MAG_BEST parameter derived by SExtrac- 
tor. The little variations of the HST/ACS PSF compared to 
the size of our targets of interest guarantees homogeneous 
photometric dataset. Table [2] lists the results from this anal- 
ysis for the 59 galaxies confirmed to be cluster members, 
using the published data listed above or the results from our 
spectroscopic survey (see below). For each target it includes 
the SExtractor ID, the coordinates (RA and DEC), and the 
photometry in each band plus errors. Due to the change of 
orientation and slightly difi^erent pointing of the individual 
HST/ACS frames, the F850LP-band covers a slightly larger 
area than the images in the other bands, and some of the 
objects were only detected in this band. 



3.2 Morphological analysis 

One of the goals of the present study is to determine the 
stellar populations of galaxies in the cluster as a function of 
their morphologies, in order to understand their origin and 
evolution. To achive this goal we performed a morphological 
analysis of the HST/ACS F850LP images, the dee pest of our 
dataset, following the prescriptions presented in ICoe et aP 
l|2006l ). where similar images were analyzed. This analysis 
recovers, for each galaxy, the Sersic index, the integrated 
magnitude, the effective radius, the ellipticity, the position 
angle, and the asymmetry index. 

For undisturbed galaxies the Sersic i ndex n correlates 
fairly well with morphological type (e.g. lAndredakis et al] 



1 19951 : [ Bell et al]l2004l ). in the sense that galaxies with high 
n have higher central concentrations (i.e., bulge dominated 
galaxies) than galaxies with low n (i.e., disk dominated 
galaxies). Less well behaved galaxies, including mergers and 
irregulars, generally do not have well defined Sersic indices. 
Fortunatelly these galaxies can generally be weeded out 
(or selected for) by m easuring their large asymmetries (e.g. 
IConselice eralll2003l ). 

In summary the morphological anal ysis consists of the 
following steps: First we use SExtractor (|Bertin et al.lll996h 
to detect galaxies in the combined F850LP-band image. 
For each detected galaxy we recover the physical and im- 
age coordinates (RA,DEC, X and Y), the integrated mag- 
nitude (mag_best), the scale-length (rso), the position an- 
gle and the ellipticity. Then, we create a postage stamp 
5xr5o on a side, for the brighest galaxy. Within that postage 
stamp, neighbouring galaxies are masked out using ellipses, 
each ellipse given a minor axis length b = 2 x rso for 
that galaxy. Using galfit ijPeng et al.l l20o3 ). the bright- 
est galaxy is fitted to a single component Sersic model 
E(i?) = Ee X exp(-K„[(7?/i?e)^/" - 1]), where is the 
effective or half-light radius, Ee is the effective surface den- 
sity, E (R) is the surface density as a function of radius and 
K — K (n) is a normalization constant. The fit is constrained 
to 0.2< n <8 and 0.3< Re <500 pixels, and the centroid 
is confined to within 2 pixels of the position derived by 
SExtractor. As initial guesses for the galfit parameters, 
we use the SExtractor output parameters. Lacking esti- 
mates for the Sersic index from SExtractor, we start all 
fits with n = 1.5. galfit convolves the model with the re- 
quired PSF. Prior to the fit, we created a PSF image by 
combining postage stamps of unsaturated point-like objects 
extracted from the combined F850LP-band image. 

Having been calculated for the brightest galaxy, the 
Sersic model is substracted from the F850LP-band im- 
age. This subtraction benefits the subsequent modelling of 
fainter nearby galaxies. We then proceed to model the sec- 
ond brightest galaxy, and continue in order of decreasing 
brightness, modelling and subtracting every galaxy in the 
SExtractor catalogue, galfit is able of simultaneous fit 
several targets, producing reasonable results, in most of 
the cases. However, in this case the fitting process is dom- 
inated by the brighter targets, and flux of the fainter ones 
is more poorly recovered. We performed extensive tests be- 
fore adopting the described procedure. In addition, we de- 
termine the asymm e try in dex (A), following the prescrip- 
tions in ICoe et al.l l|2006h . Of the 802 objects detected 
by SExtractor in the F850LP-band image brighter than 
z <26.5 mag, we applied galfit to the 241 galaxies brighter 
than z <22.5 mags. The error in the derived Sersic index is 
restricted to a range of ±0. 25, based on the simulations pre- 
sented in ICoe et al] (|2006l ) , and therefore a morphological 
segregation based on that parameter is reliable. The remain- 
ing galaxies are far beyond the depth of our spectroscopic 
data, and therefore excluded. Table |3] lists the resulting fit- 
ted parameters and their uncertainties for the 59 galaxies 
confirmed to be cluster members, using the published data 
listed above or the results from our spectroscopic survey 
(see below). It includes, for each galaxy, the SExtractor 
identification, the /'^ of the fitting procedure, the 2-band 
magnitude, the effective radius (Re), the axis ratio (a/b), the 
position angle (9), the Sersic index (us) and the Asymmetry 
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index (A) . The listed errors were esti mated based on t he re- 
sults of the simulations presented in lCoe et alj (|2006l ). The 
error in the Asymmetry index is dificult to quantify for in- 
dividual objects, however, based on the quoted simulations, 
we estimate it to be lower than ~0.1 for the galaxies of our 
sample. Of all these galaxies, only one has a strong Asymme- 
try index (~0.8). This galaxy is covered by the field-of-view 
of our IPS data, and it will be discussed later in detail. 

3.3 The Sersic indices and the Hubble sequence 

Our morphological classification scheme is based on the 
Sersic index of the galaxy profile. We consider this auto- 
matic classification scheme more straight forward, easy to 
reproduce, and less subjetive tha n the classical by-eye clas- 
sification fread lWolf et al.ll2005l . for an illustration of the 
problem). However, most of the previous studies on the mor- 
phologies of galaxies in clusters are base d on by-eye c lassi- 
fications following the Hubble scheme (jHubbld Il926h . re- 
vised by Ide Vaucouleurs et al.l (|l976l ) . Normally, they sub- 
divide the galaxies in their different Hubble types, or at least 
in three groups: pure sp heroidals (E-type), lenticulars (SO), 
and spirals (L-type) (e-g . ^Fasano et al.ll200ol : iPoggianti et al.l 
l200ll : ]ziegler et al.ll200ll : ISmail et aLh oOl'V In order to com- 
pare with previous results it is necesary to determine to 
which type our galaxies, based on their Sersic indices, cor- 
respond. 

As we indicated above, lAndredakis et aP l|l995h show 
that there is a correlation between the Sersic index and 
morphological type. However, here we need to determine 
how these p arameters correlate for our dataset. For doing 
so, we follow iFasano et al.l l|2000l ) and create a set of simu- 
lated galaxies. Their magnitudes, scale-lengths, position an- 
gles and ellipticities are extracted from the previously de- 
scribed morphological catalogue of 241 galaxies. For each 
entry in the morphological catalogue 20 galaxies were simu- 
lated, each one with a different bulge-to-disk ratio, ranging 
randomly from 0.01 (almost a pure disk) to 1000 (almost 
a pure spheroidal). Each galaxy image was created using 
galf it, with a two components Sersic model, one with ns=4 
for the bulge and another with ns = l for the disk, convolving 
them with the same PSF used in the morphological analy- 
sis. Finall y, noise was adde d foll owing the prescription pre- 
sented in ICoe et all (|2006h and ISanchez et"all (|2004cl ') for 
simulating ACS images, to resemble as much as possible the 
original F850LP-band image. 

The same morphological analysis described in the pre- 
vious section was performed for each of the simulated 
galaxy images, fitting them to a single component Sersic 
model. Finally we end up with ^5000 simulations com- 
prising the recovered Sersic index f o r a c ertain bulge-to- 
disk ratio. ISimien fc de Vaucouleur"3 (jl986l ) show the typ- 
ical bulge-to-disk ratios for the different families of galaxies 
in the revised Hubble sequence. In particular, they present 
the mean fractional luminosity of the spheroidal compo- 
nent of the galaxies and its corresponding magnitude dif- 
ference for each galaxy type. Figure |4] shows the Sersic in- 
dices of the simulated galaxies as a function of these magni- 
tude differences (Amag). The differences for various fam- 
ilies of galaxies along the Hubble sequence ar e indicated 
(extracted from ISimien &: de Vaucouleurj|l986h . Based on 
this figure we adopted the following scheme to morpholog- 



ically classify galaxies based on the Sersic indices: galaxies 
with Us >2.5 are considered spheroidals (E-type), galax- 
ies with 1.75< Us <2.5 are considered intermediate type 
(mostly SO's, but strongly contaminated by other families of 
bulge-dominated spirals, like Sa and Sb), and galaxies with 
Us <1.75 are considered spirals (L-type). 

This classification scheme is based on the fair corre- 
lation between the bulge -to-total luminosity ratio a nd th e 
bubble type of a galaxy ISimien &: de Vaucouleurs! l|l986t ). 
However we may admit that some degree of contamination is 
expected among the different morphological subgroups due 
to the dispersion in this correlation. Although we can clearly 
distiguish between pure spheroidals (E-type) and spirals (L- 
type), the subgroup of intermediate type is most probably 
less accurately defined, with considerable contaminations 
from both extreme subgroups. Another source of criticism 
about the adopted scheme is the fact that som e dwarph ellip- 
tical galaxies shown low Sersic indices (e.g.. iBarazza et al.l 
l2003h . Although in galaxy surveys without morphological 
pre-selections mos t of the galaxies with l ow Sersic indices 
are spiral galaxies lAndredakis et aP (|l995l l , the contamina- 
tion from dwarph ellipticals could be a potential problem 
in the study of galaxy clusters. To address this problem we 
perform a visual inspection of the 29 galaxies with Us <1.75 
listed in Table [2] and [3] (ie., the ones classified as spirals) us- 
ing the previous described HST/ACS images. Only 2 of these 
29 galaxies (~7%) could be classified as dwarph spheroidals 
based on this visual inspection. However both of them are 
nearby brighter galaxies, and their low surface brightness 
areas (like disk structures) are not clearly distinguished. In 
any case, the possible contamination by dwarph spheroidal 
galaxies is expected to be marginal. The remaining galaxies 
show clear evidence of disk structures. Taking into account 
these caveats we consider that the addopted classification 
scheme is valid within the context of this study. 

3.4 Spectra extraction 

To deblend and extract the integrated spectra of each indi- 
vidual galaxy in the field we used a technique for 3D crowded 
field spectroscopy developed by ourselves (ISanchez et al.l 
l2004bl : lGarcfa-Lorenzo et al ]|2005l : ISanch"e7 et al .II2OO6II. The 
technique is an extension to IFS of galf it (|Peng et al.l 
l2002h . that we named galfitSd. It consists of a deblend- 
ing of the spectrum of each object in the datacube by fit- 
ting of analytical models. IFS data may be understood as 
a set of adjacents narrow-band images, with the width of 
the spectral pixel. For each narrow-band image it is possi- 
ble to apply modelling techniques developed for 2D imaging, 
like galf it, and extract the morphological and flux infor- 
mation for each object in the field at each wavelength. The 
spectra of all the objects is extracted after repeating the 
technique through the datacube. We have already shown 
that the use of additional information to constrain the mor- 
phologic al parameters i ncreases the quality of the recovered 
spectra ISanchez et all [2OO6). With that purpose we used 
the morphological parameters previously derived from the 
HST/ACS F850LP-band image. 

We create a 2D image by coadding the flux through 
the usefull wavelength range (4650-8000 A) covered by our 
datacube. After this, we cross-checked by eye this image, 
looking for those galaxies detected in the F850LP-band for 
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Table 3. Summary of the morphological analysis. 
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Figure 5. Grayscale representation of the 2D image obtained by 
coadding the flux from the datacube through all the coveraged 
wavelength range (4650-8000 A). The circles indicate the location 
of the detected galaxies in the datacube, together with a number 
indicating the label code used to identify them. 



which we had performed a morphological analysis. We de- 
tect 41 objects (40 galaxies and 1 star) in the collapsed 2D 
image, out to 2 < 23 mag. The gravitational arc-lenses were 
excluded since their study is out the scope of the current 
work. Figure [5] shows a greyscale of the coadded 2D image 
derived from the datacube, together with an identification 
label at the location of each of the detected objects. 

As we mentioned before, we use the previously derived 
morphological parameters for each galaxy to extract its inte- 
grated spectrum by modelling each galaxy in the datacube 
with galfitSd. We use the same model that was used to 
fit the F850LP-band image, with all the morphological pa- 
rameters fixed, and fitting only the flux at each wavelength. 
For each wavelength, we created a PSF using postage stamp 
images of the only field-star in the field-of-view of our IFS 
data. As in the previous morphological analysis, the 3D fit 
was performed in a sequential way. We first extracted the 
spectrum of the brightest galaxy, the central cD, masking 
the neighbouring galaxies. Once extracted, its 3D model 
was subtracted to the original datacube, and the spectrum 
of the second brightest galaxy was extracted. The process 
continued iteratively until the faintest galaxy spectrum was 
extracted. The individual spectra are shown in Appendix 
\K\ All of them are background limited. Thus, the signal- 
to-noise ratio can be estimated from the background-noise 
level. We used the variations in the sky spectrum in each 
sky-fiber to determine the level of background noise. In this 
way, we estimated that the signal-to-noise of the extracted 
spectra range between ~1000 and ~5 per spectral pixel for 
the brightest and faintest targets. The signal-to-noise ratio 
is affected strongly by the sky-emission lines (see Fig. [T|, 
where it drops to less than ~3(t even for the brightest tar- 
gets, respectively. The effect is more severe in the wavelength 



range A > 7200 A where many sky-emission lines happen to 
be blended. 

The procedure relies on the assumption that the mor- 
phological parameters do not change substantially as a func- 
tion of the wavelength, something which it is not completely 
true, since color gradients are expected in all the galaxies. 
But, if they change, we expect this effect to have a minim 
impact in the photometry derived with galfit. Thus, if 
we use the morphological parameters derived by applying 
galfit to the F850LP-band and we force the fitting proce- 
dure to use them for deriving the photometry at any other 
band, the recovered fiux should be basically the same as 
when we let galfit fit these parameters freely. Although 
our previous experiments s hown that this is the case (e.g. 
iGarcfa-Lorenzo et al.l l2005h , we have also tested this hy- 
pothesis here, as described in Appendix[B] In any case, the 
method is not strongly s e nsitiv e to this assumption, as de- 
scribed in ISanchez et al] l|2006l ). In principle, it is possible 
to perform the extraction procedure using the parameters 
derived from the morphological analysis as an initial guess, 
and letting galfit to derive the correct ones for each wave- 
length. In this case, the recovered spectra are noiser and in 
the residual datacube there is clear evidence of not correctly 
extracted structure. 



3.5 Redshift determination 

We determined the redshift of each detected galaxy in the 
field-of-view of our IFS data by comparing its spectrum with 
syntethic models fo r single stellar populations obtained from 
iBruzual fc CharlotI l|2003l ). We will describe later in detail 
which models were used. We manually tuned the redshift 
of the model matching the different absorption features, in- 
cluding the Balmer lines (H/3, H7, H(5), the different Fe and 
Mg lines covered by our spectral range, as well as spectral 
features, like the 4000 A break. Due to the spectral resolu- 
tion of our IFS data, and the signal-to-noise of our targets, 
we do not expect to have an accuracy better than 0.001 
in the redshift determination. We derived realiable redshifts 
for 31 galaxies, 28 of them at the redshift of the cluster 
(z — 0.17±0.2), one is a foreground galaxy (z = 0.104), and 
other two are backgroud galaxies (both at 2 ~ 0.426). Table 
[4] lists the derived redshifts together with the identifications 
shown in Fig. [5] (that we will use throughout this paper) and 
the SExtractor identification listed in Table [3] to allow for 
direct comparison. 

Eleven of the 28 galaxies were studied before by 
iLe Borgne et all l| 19921 ). The redshifts derived from our data 
are coincident with the reported in this article within the 
expected errors (ctz ~0.001). We derive new redshifts for 
17 galaxies, all of them compatible with the redshift of the 
cluster. These new galaxies have been included in the Table 
[31 In total, the final sample of cluster members (28) is com- 
parable with that one of previous similar studies, but with 
the advantage that no sample has been selected: our sample 
is complete within the field-of-view of our IFS data down 
to our detection limit. Of these galaxies, only 6 have been 
previous studied with enough quality spectr oscopy to de- 
rive t he properties of their stellar populations l|Ziegler et al.l 
120011 ). Therefore, our sample will increase significantly our 
knowledge of the core of the cluster. 
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Table 4. Summary of the parameters derived from the spectra 
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Uncertain values are written in italic. 



3.6 The color-magnitude diagram 



3.6.1 Morphological segregation 



Figure [S] shows the rest-frame B-V vs. Mv^ color-magnitude 
diagram for those cluster members listed in Table O de- 
rived from the HST/ACS based photometric data described 
above. The absolute magnitudes were corrected for cosmo- 
logical effects and we applied a fc-correction based on the 
same SSP synthetic model. For each synthetic SSP model 
(described below) we derived the expected B-V colors at 
the average redshift of the cluster and in the rest-frame, de- 
riving at the same time the V-band fc-correction. An aver- 
age relation between the observed and rest-frame colors and 
the fc-correction was derived, and applied to the data shown 
in Fig. [6| Red solid circles indicate galaxies with Sersic in- 
dex n >2.5 (E-type), while green stars indicate those with 
n <1.75 (L-type). Orange solid squares indicate galaxies 
with Sersic indices in between those values, corresponding 
to intermediate galaxies. The cD galaxy is indicated with 
an orange open cross, being the brightest object in the field. 
We prefer to label it with a different symbol since it is well 
known that cD galaxies do not follow a Sersic profile of any 
index, showing a flatt er and more e xtended hght profile in 
the outer regions fe.g. IOemleilll976l) . Open circles indicate 
the galaxies within the field-of-view of our IFS, with red- 
shifts listed in Table (4] The big open square indicates the 
only galaxy with a high Asymmetry index. 



To our knowledge this is the first time that a morphologi- 
cally segregated color-magnitude diagram is shown for con- 
firmed cluster members. In previous works (e.g. Smail et al. 
2001) only cluster galaxy candidates have been shown. The 
morphologies, based on the Sersic indices listed in Table |3l 
break down as: 11 E-type, 19 intermediate and 29 L-type. 
If we consider as E-type galaxies to be all the objects with 
n >1.75, they represent only a ~50% of the total. For the 
complete subsample of galaxies covered within the IFU field- 
of-view, the morphologies breaks down as: 7 E-type, 8 inter- 
mediate and 13 L-type. Contrary to previous results (e.g. 
IZiegler et aLll200ll : ISmail et al. 2001 1 this indicates that the 
cluster population is not dominated by early-type galaxies, 
and there is an almost parity between early- and late-type 
galaxies. 

Our approach to derive the morphologies is substan- 
tially different to in previous attempts, based on visual 
classifications. There is a considerable debate over the re- 
liability of visual classifications of early-type galaxies (e.g . 
iDressler et al. I ll997l : lAndreon|[l99i : iFabricant etlo] |200G| ) . 
and the ability to compare them with classifications based 
on profile analysis, like concentration and/or Sersic indices. 
However, since there are no previ ous similar stu d ies, w e must 
compare with the existing ones. IZiegler et all |20o3) show 
visual morphological classifications of 19 of their 48 cluster 
members, based on WFPC/HST imaging. Their classifica- 
tion was performed by W. Couch, who classified them as 8 E, 
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1 E/SO, 5 so, 3 SB/0, 1 Sa and 1 Sab. Taking together all the 
E-type galaxies (E, E/SO and SO) , they account for a ~73% 
of the total sample. ISmail et al.l ()200lh presented morpho- 
logical classification for 81 cluster candidates (33 with con- 
firmed redshifts), based on the same analysis as lZiegler et al.l 
l|200ir i. Their morphologies break down as: 31 E, 3 E/SO, 28 
SO, 5 Sa and 14 late-type spirals. Taking together all the E- 
type galaxies, they account for a --^76% of the total sample. 
Even taking into account the different methods used to de- 
rive the morphologies , and the diffi c ulties to co mpare them, 
it is c lear that both IZiegler et al.l (1200 ll ) and ISmail et al.l 
l|200ll ) found fewer L-type galaxies than what we are find- 
ing.^ 

IZiegler et al.l (|200ll ) selected galaxie s using the pre- 
viousl y existing redshift catalogue from iLe Borgne et al.l 
and defining a region in the UBVI colour space 
which was occupied by cluster members. Finally, only galax- 
ies brighter than / <18.8 mag were targeted by their spec- 
troscopic survey. Although it was claimed that this selec- 
tion places negligible restrictions on the stellar populations 
of the selected galaxies, and was made with the aim of re- 
jecting the majority of background galaxies, it may have 
had an effect on the morphological segregation. The bright- 
ness cut in the 7-band implies that all the observed galaxies 
are brighter than My < ~-19.5 mag, asuming an average 
V — I color of ~1.5 mag (|Fukugita et al.lll995l ). If we ap- 
ply such a cut in our sample (see Fig. the morpholo- 
gies break down as: 9 (7) E-type, 14 (5) SO's and 11 (2) 
L-type for the full (IFU) sample. This means a fraction of 
^68% (86 %) of E-type g alaxie s, in agreement with the re- 
sults from [Ziegler et al] (|200ll ). Therefore, although visual 
classifications are not fully reproductible and may produce 
discrepancies in the morphological classification of individ- 
ual objects, they agree with those ones based on Sersic in- 
dices on a statistical base. 

I Small et al.l (|200ll ) selected 81 cluster candidates based 
on a deep Ks-hand survey {Ks <19 mag), discarding those 
objects with clear stellar morphology, in too crowded areas 
and /or showing lensed morphologies. Assuming a z — Ks 
color of ~3 mag for the reddest objects, based on their own 
results for I — Ks an d an average z — I color of ~0.5 mag 
iFukugita et al.lll995l ). they selected galaxies brighter than 
z <~22 mag, a limiting magnitude similar to that of our 
IFU data. Therefore, a difference in the depth of both sur- 
veys is not the reason for the differences in the morphologi- 
cal segregation. We matched their cluster candidate sample 
with our full sample (listed in Table [3]) , finding 39 coinci- 
dences. For this subsample of confirmed cluster members 
the morphological segregation based on their analysis re- 
mains similar, with a ~82% of E-type galaxies. Therefore, 
possible contaminations from non cluster members are not 
the cause of the discrepancy. We examined the Sersic in- 
dices of those 32 E-type galaxies. Of them, 10 have Sersic 
indices lower than n <1.75, being L-type galaxies following 
our scheme, and only 7 have n >2.5. Thus, the discrepancy is 
most probably due to the different criteria adopted for clas- 
sifing the objects, illustrating the risks of direct comparisons 
between both methods. Statistically speaking visual classifi- 
cation gives the same distribution of galaxy types as profile 
fitting: on average, their E-type galaxies have larger Sersic 
indices (n ~2.6) than their L-type ones (n ~1.7). However, 
it is well known that one-to-one cross-checking both meth- 
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Figure 6. Rest-frame B-V color distribution as a function of 
absolute magnitude My for those cluster members listed in Table 
[2] Red solid circles indicate galaxies with Sersic index n >2.5 (E- 
type), while green stars correspond to objects with with n <1.75 
(L-type). Orange solid squares indicate lenticular galaxies, with 
Sersic indices in between both values. The cD galaxy is indicated 
with an orange open cross, being the brightest object in the field. 
Open circles indicate the 28 galaxies within the field-of-view of 
our IFS. The big open square indicates the only galaxy with a 
high Asymmetry index. 

ods give different results (e.g. IWolf et al.ll2005l ). and profile 
analysis produces a more realiable and objetive scheme for 
classifying objects. 

3.6.2 Analysis of the CM diagram 

The Color-Magnitude (CM) diagram contains extremly use- 
ful information to understand the nature of the cluster pop- 
ulation. The B — V color is related to both age (strongly) 
and metallicity (weakly) of the dominant stellar population 
of the galaxy. On the other hand, the V^-band absolute mag- 
nitude can be related to the stellar mass. Figure[6]shows a bi- 
modal distribution of early- and late-type galaxies in the CM 
diagram. E-type galaxies are (mostly) located in the brighter 
end of the distribution (no one is fainter than Mv <~ —20.5 
mag), covering a narrow range of colors (1.9< B — V <2.2). 
This is a well known feature of early type galaxies in clusters 
and it has been used fo r decades to select cluster candidate s 
(e.g. IZiegler et allbOOll : ISanchez fc Gonzalez-Serranoll2002l ) , 
and to trace the evolution of the older popu lation of cluster 
galaxies (e.g. lAragon-Salamanca et al ] | 19931 ). 

On the other hand, L-type galaxies are distributed 
across a wide range of magnitudes with the brightests al- 
most as red as the early type galaxies, and the faintests 
bluer and covering a wider range of colors. For galaxies 
brighter than Mv < —19.7 mag there are no significant 
differences between the core members (galaxies in the IFU 
field-of-view), and the rest of the galaxies of the sample. For 
galaxies fainter than Mv > —19.7, it is found that the clus- 
ter core galaxies cover a wider range of colors, with more red 
galaxies. However, this difference may be due to a selection 
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effect of the galaxy samples obtained from the literature. 
As indicated in Section 4.3 two of the 29 galaxies classified 
as L-type could be spheroidals on the basis of a visual in- 
spection of the HST/ACS images. Both of them have red 
colors, B — V ~0.9 mag. However they are not the redder 
L-type galaxies in their magnitude range. It is interesting 
to note here the location of the galaxy number 38 (IFUid), 
the only one showing a high asymmetry, being the bluest 
galaxy of our sample. This galaxy, at the East in Fig. 2 
and 3, is clearly identified as a merger in its lastest stages, 
showing a clear tail to the South- West. Its spectrum. Fig. 
lAll shows strong H/3, [OIII]4959,5007 and Ha emission lines 
and weaker [NII]6548,6584, indicating a large amount of star 
formation. Its blue colors and spectrum also indicate the 
presence of an extremely young stellar population, far away 
from the global trend of the rest of the cluster members, 
as we will discuss below. Contrary to what is expected this 
galaxy is just at the nominal redshift of the cluster {z = 0.17, 
based on the emission lines), and considering its projected 
distance to the cluster center (assumed to be the cD), it 
is not likely to be a fly-by object. To our knowledge it has 
escaped the attemption of previous studies, most probably 
because its morphology and color automatically exclude it 
as a cluster candidate. Our IFU data lack the required spa- 
tial and spectral resolution to perform a more detail study 
of its properties. 

A brief comparison of the color-magnitude distribution 
of the different kind of galaxies in the cluster, segregated by 
their morphology, indicates that E-type galaxies are more 
meissive and have older and more metal rich stellar popu- 
lations. L-type galaxies cover a wider range of masses, ages 
and metallicities, with the more massive ones being similar 
to the E-type galaxies, and the less massive ones younger 
and metal poorer. This result indicates that a bulk star 
formation process at high redshift may not be the unique 
mechanism to explain the population of the galaxies in the 
cluster. We will investigate it further in the next sections. 

3.7 Galaxies ages and metallicities 

We obtained spectra of the individual galaxies in the cluster 
with the aim of comparing their global spectroscopic proper- 
ties with their morphological ones in a relative way, ie., only 
within the cluster. It is well known that the spectral energy 
distribution (SED) of simple stellar populations (chemically 
homogeneous and coeval stellar systems) depends on a set 
of first principles (e.g., initial mass function, star forma- 
tion rate, stellar isochrones, element abundance ratios), to 
make it possible to generate synthetic stellar populations 
from them. This technique, known as evolutionary synthesis 
modeling (e.g. Tinsley 1980), has been widely used to unveil 
the stellar population content of galaxies by reconciling the 
observed spectral energy distributions with those predicted 
by the theoretical framework. Unfortunately the variation 
of different physical quantities governing the evolution of 
stellar populations produce similar effects in the integrated 
light of those systems, leading to a situation in which the 
observational data is affected by undesirable degeneracies, 
like the widely mentioned one between age and metallicity 
(e.g. O'ConneU 1976; Aaronson et al. 1978, Worthey 1994a). 
The most widely used technique to compare galaxies with 
models is to measure certain line strength indices, on the 
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Figure 7. Example of the method used to derive the H/3 and 
Mg6 line indices. The plot shows a detail of the spectrum of the 
galaxy number 2 (IFUid), zoomed in onto the region of interest. 
For each absorption line the boxes indicate the regions were the 
adjacent continuum are estimated (blue and red boxes) together 
with the region where the absorption is measured. The orange 
line indicates the adopted slope for the continuum. The strong 
feature at 5577A is the residual of the non perfect subtraction of 
a night sky emission line. 



Lick/IDS index system (Burstein et al. 1984; Faber et al. 
1985; Burstein, Faber & Gonzalez 1986; Gorgas et al. 1993; 
Worthey et al. 1994). People generally try to use a combina- 
tion of indices that are most orthogonal in parameter space 
(e.g. age and metallicity; although see Cardiel et al. 2003 for 
a discussion of the impact of random errors in the selection 
of the best combination of indices). Following the Lick/IDS 
definitions, for our study we have chosen H/3 as a primary 
age indicator, and the Mgb index as a metallicity indicator. 
We use the Mgb index rather tha n the more convent ional 
[MgFe] one, because both we and IZiegler et ahl |20o3) are 
unable to derive realiable Fe indices for most of the galax- 
ies. This diagram has the drawback that the models based 
on solar abundance ratios do not match the Mg/Fe rati o 
observed in early-type galaxies (|Thomas et al.|[l999l . |2003 ). 
It is also well known that the ages derived from the grid 
are clearly uncertain in an absolute sense. However, we are 
primarily interested in a relative comparison of the ages and 
metallicities of the galaxies in the cluster, and in this sense 
the diagram is valid. 

As we will discuss below, a proper estimate of those 
indices requires high signal-to-noise spectra, and may be 
strongly affected by the instrumental and systemic veloc- 
ity dispersions, the presence of weak emission lines that fill 
the absorptions (specially important in the case of H/3), and 
the presence of high residuals corresponding to the removal 
of nearby sky emission lines. 



3. 7. 1 Analysis of line-strength indices 

We derive the H/3 and Mg6 equivalent widths for our IFU 
sample using the bandpass defi nitions from the Lick index 
system revised by iTrager et al. 

I 11998), shifted to the red- 
shift of each object. We used our own routines and checked 
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Figure 8. Age - metallicity diagnostic diagram using H/3 as a 
primary age indicator and Mgfe as a metallicity indicator. The 
symbols are similar to those ones in Fig|6] Overplotted is a grid 
of SSP models from Worthey (1994b). 'Horizontal' lines follow 
constant age (Gyr), 'vertical' lines constant metallicity ([M/H]). 
Each line is labeled with its corresponding age/metallicity. 

the results with index^El. Figure [7] illustrates the procedure, 
showing a detail of the spectrum of galaxy number 2 and 
the actual bandwidths used to derive each of the indices. 
The derived results are listed in Tabled together with the 
error estimates (computed as explained below) . We only list 
the positive values, since a negative equivalent width means 
emission, and not absorption. In most of the cases these 
negative values are found in faint targets, and are simply 
due to low signal-to-noise of the spectra. However, in a few 
cases they indicate real emission lines, like the strong emis- 
sion lines of galaxy number 38 (the merger) or fainter ones, 
like in galaxies numbers 18 and 22. In these three cases it is 
possible to derive realiable Mg6 indices, but unreliable H/3 
ones. On the other hand, there are a few objects (16, 20, 
21 and 25), with too low Mgfc equivalent widths for their 
H/3 ones. A visual inspection of their spectra, and a trace 
back to the raw data shows that they show cross-talk con- 
tamination from ThAr emission lines obtained in adjacent 
calibration fibers, which unfortunally lies in the regime of 
the Mgfo absorption. For other four galaxies (number 1, 7, 
8 and 15) both H/? and Mgb seem to be too low. Galaxy 
number 1 is the cD, a kind of galaxy with wel l known devi- 
ation s of the metacillity sensitive indices (e.g. ICardiel et al.l 
1 19981 ). A visual inspection of its spectrum in the H/3 region 
indicates possible emission that fill the absorption. A similar 
inspection of the spectra of galaxies number 7 and 8 shows 
clear contaminations of sky emission line residuals in the 
adjacent regions to the H/3 absorption band. Finally, galaxy 
number 15 has an uncertain redshift determination, which 
strongly affects the derived indices. All these objects have 
been excluded from further analysi s based on line indi ces. 
Of the 6 galaxies in common with IZiegler et al.l (|200ll ) , 5 
have line-strength indices in common. In both cases, the 
strengths of H/3 and Mgfe agree within ±0.6A, which is far 
beyond the estimated errors of both studies. However such 
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differences are expected on the basis of the different meth- 
ods used to derive the strength of the line indices that may 
strongly affect the results. 

We estimated the errors using Monte Carlo simulations 
of each spectrum, in the following way: first we derive for 
each spectrum the best matched single stellar population 
model, as we will describe in the next section. This model is 
included in Fig. lAll Then, we substract the model from the 
spectrum, and derive the absolute value of the difference for 
each wavelength. Finally, we smooth this residual spectrum 
with a 10 pixels width median-box filtering. This smooth- 
ing was performed to reproduce as much as possible the 
noise cross-correlation pattern between adjacent pixels. The 
width was selected to match the spectral resolution, since 
we consider that this is the range of pixels where the cross- 
correlation of the noise is stronger. We consider that this 
rather complex method to estimate the errors is a better 
representation of the variance than the much simplier ap- 
proach that would be to create a variance spectrum based 
on the estimated accuracy of our spectrophotometry (Sec. 
2.1 and Appendix B). The reason for that is that using this 
latter approach we would not take local effects in the vari- 
ance, like unperfections in the sky subtraction, cross-talk 
contaminations and unaccuracies on the reduction process. 
Using the adopted variance and our original spectrum we 
created 100 simulated spectra. The standard-deviation of 
the resulting equivalent widths are assumed to be the errors 
of our measurements. Line indices, for galaxies with velocity 
dispersions larger than 200 Km/s (derived from the fitting 
to SSP models, as described below), are finally corrected as 
in Jablonka, Gorgas and Goudfrooij (2002) by using similar 
coefficients (Jablonka, private communication). 

The final useful dataset comprises line indices for 26 
objects, 12 of them corresponding to our IFU sample. Fig- 
ure [8] shows the Age - metallicity diagnostic diagram using 
H/3 as a primary age indicator and Mgfe as a metallicity 
indicator. The symbols are similar to those in Fig|BJ includ- 
ing both our determined equ ivalent widths and those ob- 
tained by'Zioglcr et al.' ('2001'). Overplotted is a grid of SSP 
models from Worthey (1994b). 'Horizontal' lines follow con- 
stant age, 'vertical' lines constant metallicity. Even though 
there are in the literature more recent evolutionary synthe- 
sis model predictions suitable for our purposes (e.g. Bruzual 
& Chariot 2003, Vazdekis et al. 2003, Thomas, Maraston & 
Kron 2004), here we are using the classical Worthey (1994b) 
models because we are just interested in a differential anal- 
ysis of our data (i.e. the absolute derived ages and metal- 
licities should not be taken literally) and we can perform a 
direct co mparison with the re sults published by Ziegler et 
al. (2001). IZiegler et al] l|200ll ) already noticed that E-type 
galaxies in this cluster cover a wide range of parameters 
within Fig [51 mu ch wider than similar samples in the lo- 
cal universe (e.g. iKuntschner et al]|200ll ). We confirm this 
conclusion, since the new galaxies have increased the range 
of covered parameters. The bulk of the galaxies show lower 
metallicities and younger ages at a fixed metallicity com- 
pared to the ones in the local universe. They suggested that 
these differences may be due to variations in the average 
masses of the two galaxy samples, in the sense that more 
massive galaxies tend to be older. 

The IZiegler et al.l l|200ll ) sample lacks enough L-type 
galaxies to draw a clear conclusion on the dependence of 
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these parameters on the morphology. With the new sample 
it is possible to see that L-type galaxies are found towards 
younger ages, on average, than E-type galaxies at any fixed 
metallicity. On the other hand they seem to cover the same 
range of metallicities at any fixed age. These conclusions 
have to be taken with care, since most of our L-type galax- 
ies are, at the same time, less luminous and therefore less 
massive than our E-type ones (Fig. [6]). Therefore, if there is 
a correlation betwe en the age and the mass, as suggested by 
IZiegler et all |20o3), this trend may affect both the E-type 
and the L-type galaxies, and explain our results. We will 
further investigate this possibility below. 



3.7.2 Fitting to Single Stellar Population models 

ICardiel et al.l (|2003t ) already showed how sensitive the 
derivation of physical parameters from line indices is to the 
signal-to-noise of the spectra. They concluded that using the 
full spectrum, when possible, one can derive better param- 
eters, for a similar signal-to-noise. In the previous section 
we have seen that even for high signal-to-noise spectra local 
contaminations from emission lines, both from the object 
and/or residuals from the sky or cross-talk from calibration 
fibers can strongly affect the derived indices. 

To extract the maximum information for the maximum 
number of objects we fitted each spectrum from our IFU 
sample with synthetic SSP models. Models were created us- 
ii ig the GISSEL code (jBruzual fc Charlotll2003h. a ssuming 
a lSalpeteil l| 19551 ) IMF fthe use of a lChabrieJ (|2003l ') will not 
modify the results), for different ages and metallicities. We 
create 72 models covering a discrete grid of 12 ages (5 Myr, 
25 Myr, 100 Myr, 290 Myr, 640 Myr, 0.9 Gyr, 1.4 Gyr, 2.5 
Gyr, 5 Gyr, 11 Gyr, 13 Gyr, and 17 Gyr), and 6 metallicities 
{Z =0.0001, 0.0004, 0.004, 0.008, 0.02 and 0.05). Each spec- 
trum was then fitted t o each of the 72 models using FIT3D 
ijSanchez et al.l l2006bl ). This code resamples the model to 
the resolution of the data, convolve it with a certain veloc- 
ity dispersion, shift it to the redshift of the target and scale 
it to match the dataset by a x^/v minimization scheme. 
Therefore, for each SSP model, these are the three free pa- 
rameters in the fitting process: velocity dispersion, redshift 
and a scaling factor related to the luminosity. The redshift 
of the object was fixed to reduce the degrees of freedom. 
Only the spectral region at wavelengths bluer than 7200 A 
was used in this analysis. At redder wavelengths the strong 
night sky emission-lines and the telluric absorptions have 
not properly been corrected for in our reduction process, and 
strong residuals are seen in the spectra. Wavelength regions 
possibly affected by imperfect sky emission line subtraction 
were also masked. We derived the velocity dispersions for 
each object spectrum in the fitting process, including both 
the instrinsic and instrumental ones. We cannot perform an 
accurate correction of the instrumental effects cause we lack 
of the proper comparison observations of velocity calibra- 
tion stars during the night. This prevents us to use them as 
a mass indicator (e.g. IZiegler et al.ll200ll ). using them only 
to correct the line-indices strength, when required. The cor- 
rected valocity dispersions range between ~100-400 km s~^, 
and despite the lack of accuracy of our derive d velocity dis- 
persio ns they agree with the ones derived by IZiegler et al.l 
lj200ll ). for the 6 objects in common. 

We obtain for each spectrum a grid of 72 /v, each 
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Figure 9. Greyscale image of the 2D distribution of /v derived 
from the fitting of a single galaxy spectrum to a grid of SSP 
models of different ages and metaUicities. Darls colors indicate 
lower 'x^ /y (ie., better fits), while light colors indicate higher vales 
(ie., worse fits). The of the four best fits are labeled. The 

plotted example corresponds to the second brightest galaxy in our 
field-of-view (id=2 in Table Hll 



one corresponding to an age- metallicity pair model. Figure 
[9] illustrates the result of the fitting technique. It shows a 
greyscale of the /v map derived after fitting the spec- 
trum of galaxy number 2 to the grid of models, each one 
corresponding to an age- metallicity pair. Similar results are 
found for the rest of the objects. In the fitting procedure 
we assumed a simple root-mean-square of the flux as an es- 
timation of the uncertainties in the spectra, to derive the 
X^ /v. Therefore these values are not valid to compare the 
results of the fit for different objects, and just to perform 
internal comparisons between the results of the fit for the 
same object. The X^ /v of the four best fits are indicated 
in the figure. In general it is not possible to distinguish 
which is the best fitting model among the best four ones. 
Therefore, we estimate the age and metallicity of each ob- 
ject by averaging the ages and metallicities of the four best 
models, weighted by the inverse of the corresponding x^ /v 
values. Errors for the estimated parameters are derived by 
Monte Carlo simulations: for each spectrum we estimate the 
variance spectrum as described above. Then we created 30 
simulated spectra, using as input the original spectrum and 
adding noise following the estimated variance. For each of 
the 30 simulated spectra we run the fitting process over the 
72 models (i.e., 2160 fits per input spectrum), and the age 
and metallicity are derived following the recipe described 
above. The standard deviation of the recovered distribution 
of ages and metallicities for the simulated spectra is consid- 
ered as a good estimate of the error of the measured param- 
eters for the input spectrum. The final errors for the ages 
and metallicities are computed coadding quadratically these 
estimated errors to the standard deviation of the values for 
the four best models. The results from this procedure are 
listed in Tabled 

The method relies on certain assumptions: (i) we con- 
sider that the use of the full spectrum is a better method 
to recover the age and the metallicity of the main stellar 
population than just using the line indices and (ii) we con- 
sider that our spectrophotometric calibration is stable from 
blue-to-red, i.e., that the shape of the spectra is well recov- 
ered. This is a sensitive assumption, since the whole method 
relies on it. To test its validity we compared our integrated 
spectra with the available broad-band photometry (Fig. lAl| ). 
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The photometric dataset, derived from the HST/ACS data, 
covers the full optical wavelength range, from '^4700A to 
~9000A, and in particular the wavelength range of our spec- 
troscopic data. In most of the cases our spectroscopic data 
match pretty well with the broad-band photometric ones. 
Even more, the extrapolation of the data using the best fit- 
ted SSP synthetic model reproduces also the expected pho- 
tometry in wavelength ranges not covered by our spectro- 
scopic data in most of the cases. Only galaxies number 9, 
11, 16, 18, 22 and 24 show a slight missmatch (~2(t) be- 
tween the extrapolated flux and the zsso-band photometry. 
If we take into account that the spectrophotometry seems 
to match for the rest of the bands and for the rest of the 
objects, we have strong arguments to consider that is also 
valid for thes ones: We used an IFU, and therefore the flux 
calibration was performed only once for the full dataset. 

Finally we assume that (iii) the adopted grid of mod- 
els comprises a good representation of the expected stellar 
populations. We tried not to pre-judge the composition of 
the stellar population and cover the wider possible range of 
ages and metallicities. On the other hand, we must limit the 
number of input models for a reasonable performance of the 
fitting procedure. This imposes limitations on the sampling 
of the grid, and therefore, in the accuracy of the recovered 
age and metallicities. We must take that into account when 
drawing conclusions from our modelling. 
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Figure 10. Distribution of the metallicities as a function of age 
derived by the fitting procedure for the 28 galaxies our IFS sam- 
ple. The symbols are similar to those ones in Fig|6] The size of 
the symbols indicate the mass, with larger symbols for more mas- 
sive objects. All the E-type galaxies (n >2.5), and most of the 
n > 1.75 ones have stellar populations older than ~9 Gyr and a 
wide range of metallicities, 0.02 < Z < 0.05. 



3.7.3 The Age-Metallicity Distribution 

Figure [TD] shows the distribution of the metallicities as a 
function of age derived using the previously described fitting 
procedure for the galaxies in our IFU sample. The symbols 
are similar to the ones in Fig|S] The size of the symbols in- 
dicates the mass (see below for a description of how it was 
estimated), with larger symbols for more massive objects. 
The merging galaxy (IFUid =38), is not visible in this plot, 
since the derived age for its stellar population, ~250Myr, is 
far from the average of the rest of the sample. This result is 
expected due its B — V color, and the strong star formation 
process that is derived from the strength of the emission 
lines on its spectrum. In this particular case a SSP model 
is not a good enough to describe the expected mix of pop- 
ulations resulting from a merging process, with at least an 
underlying stellar population, and a younger one resulting 
from the induced star formation process. 

If we compare this flgure with the H/3-Mgb distribution 
shown in Fig. [8] evident differences can be seen. However, a 
direct comparison has to be done with care. First, the ab- 
solute ages and metallicities derived from both figures can- 
not be compared, since the y come from com pletely different 
SSP models. We used the IWorthevI l|l994bD models in the 
previous figure to be consistent with Ziegler et al.l (|200ll ). 
while the models used in the fitting procedure are derived 
from the GISSEL code (|Bruzual fc Cliarlot|[2003 l. However, 
relative comparisons between the different galaxies are pos- 
sible, in principle. One needs to take into account that in 
Fig[S]we only represent 12 of the 27 objects shown in the 
current plot, i.e., most ly the brightest one s. On the other 
hand, the objects from IZiegler et al.l (|200ll ) are not plotted 
in the current figure. However, the conclusions about the 
stellar populations of individual galaxies derived from both 
figures agree qualitatively: the galaxies with young and old 



stellar populations have young and old stellar populations 
in both figures, and the galaxies richer and poorer in metals 
are richer and poorer in both figures too. 

We consider the current figure a more accurate repre- 
sentation of the real distribution of ages and metallicities 
of the galaxies in the cluster (or at least in the core). The 
plotted parameters were derived for the IFU sample, which 
is complete to our detection limit, without any pre-selection 
and any a posteriori exclusion. In general it is found that the 
galaxies with older stellar populations covers a wider range 
in metallicities, being more metal rich, on average, than the 
galaxies with younger stellar populations. This result does 
not agree with the trend between these two parameter found 
in field ellipticals, where galaxies with younger stell ar pop- 
ulations are more metal rich (e.g. iTrager et aLlbOOd ). How- 
ever, our sample comprises a large fraction of L-type galax- 
ies. As in Fig [5] we found a clear different behaviour for the 
E-type and L-type galaxies. All the E-type galaxies with 
n >2.5 (spheroidals) , and most of the n > 1.75 ones (in- 
termediate galaxies) have stellar populations older than ^9 
Gyr, covering a wide range of metallicities, 0.02 < Z < 0.05. 
On the other hand, the L-type galaxies (n < 1.75), cover a 
much wider range of ages, from ~3 Gyr to ~13 Gyr (keep- 
ing in mind that the absolute ages have no meaning), and 
metallicities, 0.001 < Z < 0.05. Interestingly, the L-type 
galaxies with the largest metallicities {Z >0.015) seems to 
have similar ages as the E-type ones, and the L-type galaxies 
with smaller metallicities cover a wider range of ages. The 
most metal rich galaxies are the more massive ones (with 
larger symbols), and therefore, this result may indicate that 
more massive galaxies tend to be older whatever their mor- 
phologies. We will discuss that in detail below. This result 
is certainly not clear from Fig|8l 
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Figure 11. Distribution of the stellar masses for the galaxies of 
our sample derived from the results of the fitting procedure to 
SSP models along the stellar masses derived using the Bell & de 
Jong (2001) approximation. 
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3.8 Stellar populations and Masses of the 
Galaxies 

In previous sections we have already seen that the masses of 
the galaxies are directly related to the nature of their stellar 
populations. The Color-Magnitude diagram, Fig. |6l shows 
that less luminous (and massive) galaxies show a wider range 
of colors than more luminous objects. The comparison of the 
H/3-Mg& distribution of this cluster with that of the local uni- 
verse, including more massive galaxies, makes [Ziegler et al.l 
l|200ll ) think that there might be a trend of metallicities/ages 
with the mass of the objects. Figure [10] shows that more 
massive galaxies are, on average, older and more metal rich. 
We investigate here the possible relations between the stel- 
lar populations and the morphology with the mass of the 
galaxies. 

It is well known that the masses of spheroidal or early- 
type galaxies are directly related to their stellar velocity dis- 
persions (cr), due to their nature. This parameter, directly 
accesible from the spectra, has indeed been used as a mass 
indicator for decades, leading to a set of well establish scaling 
relations when comparing it with different properties of the 
galaxies. Among these relations the most important ones are 
the lumi nosity-velocity dispersio n relation or Faber- Jackson 
relation (|Faber fc Jacksonlll976l ). the surface brightness, ef- 
fective radi us and velocity dispe r sion relation or Fundamen- 
tal plane llDressler etM 1 19871 : iDiorgovski fc David Il987l : 



iBender et al.l 19921 ) and the Mg-yelocity dispersion relatio n 
or Mgb-a relation l|Bender et al.lll993l : ICoiless et al.lll999^ . 
These scaling relations have been used to compare between 
different p opulations of early -type galaxies at different red- 
shifts (e.g. IZiegler et al]|200ll ). to trace the evolution of this 
kind of galaxies. For late-type galaxies similar relations are 
found when using the maximum rotating velocity, a param- 
eter directly related with the mass of the rotating disk. 
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Figure 12. Distribution of ages along stellar masses for the our 
galaxy sample. The symbols are similar to those ones in Fig|6l 
There is a clear morphological segregation. L-type galaxes are 
less massive than 10^^ Mq, and have a wider range of ages than 
E-type galaxies. 



All these relations have a direct correspondence with 
relations between the mass and the different properties of 
the galaxies: The Faber- Jackson relation can be understood 
as a Luminosity-Mass relation, the Fundamental plane as 
a Mass-scale relation and the Mg-o- as a Metallicity-Mass 
relation. Our spectroscopic data lack the required resolu- 
tion and signal-to-noise for a proper determination of the 
velocity dispersion for most of the galaxies, therefore, we 
prefer to investigate these relations using the stellar masses 
of the galaxies. This approach has the advantage that we 
can explore at the same time all the galaxies of the sample 
independent of their morphology. 

We estimat e the stellar ma s ses u sing the average re- 
lation found by iBell fc de Jond (jiooj) between the M/L 
for the B-band luminosity and the B — V color. For doing 
that we derive the rest-frame B — V colors of our objects 
using the best fitted SSP model described before, synthesiz- 
ing the magnitudes by convolving the model with the filter 
responses, and obtaining the _B-band absolute magnitude 
from our photometric data. It i s important to note that the 
masses of lBell fc de Jongj (|200ll ) are just relative values, cal- 
ibrated under the assumption that spiral galaxies are in a 
maximum-disk situation. However, this will not affect our 
results since we are interested in relative comparisons be- 
tween difi'erent families of galaxies. The derived masses are 
listed in Table |4] The errors listed in the table were derived 
from the photometric errors in the B-band magnitudes listed 
in in Table [2] and they do not include any other possible 
source of error (like the unaccuracy of the fitting process or 
the derived synthetic colors) . However, we consider that this 
error dominates any other possible one. The masses of the 
objects range between ~10® and ~10^^ M©, being the cD 
galaxy the most massive galaxy and the merging galaxy the 
least massive one. 



3.8.1 The Age-Mass relation 

Figure [12] shows the age distribution as a function of stellar 
masses for the galaxies of our IFU sample. The symbols are 
similar to those in Fig[6] The division of galaxies in terms 
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Figure 13. Distribution of Sersic indices as a function of the 
stellar masses for the galaxies of our sample. The symbols are 
similar to those in Fig|6] There is a clear trend in the sense that 
E-type galaxies are more massive than L-type ones. 
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Figure 14. Distribution of the metallicities as a function of stellar 
mass for the galaxies of our sample. The symbols are similar to 
those ones in Fig[6] There is a clear correlation between both 
parameters, in the sense that more massive galaxies are more 
metal rich than less massive ones. 



of morphologies is remarkably similar to that of the Color- 
Magnitude relation (FigO, although in the current plot the 
segregation between early- and late-type galaxies is much 
stronger. Most of the E-type galaxies (n > 1.75) and all 
the ones with n > 2.5 are more massive than ~7 10^" Mq 
and have stellar populations older than ~9 Gyr. On the 
other hand, all the L-type galaxies are less massive than 
this limit, and the ages of their stellar populations spread 
towards younger ages. There is no obvious trend of age with 
mass, but there is a lack of young stellar populations in 
massive galaxies. Since all of them are E-type ones, it is 
not possible to know if this deficit is due to the mass or 
the morphology of the galaxies. However such a trend it 
is expected, a t least for field early-type galaxies. Recently 
iPeletier et al.l (|2006l ) reported a trend between the age and 
the mass for field e arly-type galaxies based on data from the 
SAURON project (|Bacon et al .11200 lI V 

Figure [13] illustrates this effect. It shows the distribu- 
tion of Sersic indices as a function of the stellar masses for 
the galaxies of our IFU sample. The symbols are similar to 
those in Fig[6l There is a strong correlation between both 
parameters, in the sence that E-type galaxies are massive 
and L-type ones less massive. The relation is clearer when 
masking the cD and the merger galaxies, for which the pro- 
files are not well represented by the adopted model. 

Similar results were already presented bv lZiegler et al.l 

(|200ll ) for the early-type galaxies of A2218 in the outer re- 
gions. They compared the H/3 index with the stellar veloc- 
ity dispersion finding that the fainter galaxies shows a larger 
scatter in their If/3 values, indicating a wide range of star for- 
mation histories. Our survey goes deeper, sampling galaxies 
much fainter (and l ess massive) that those ones studied by 
IZiegler et all (|200lh . extending their conclusions to fainter 
regimes. 



3.8.2 The Metallicity-Mass relation 

As we pointed out above, the Mg&-(7 relation can be under- 
stood as a metallicity-mass relation. Figure [2] shows the dis- 
tribution of the metallicities as a function of the stellar mass 
for the galaxies of our IFU sample. The symbols are similar 



to those ones in Fig[6l As expected there is a clear corre- 
lation between both parameters (r =0.997,P >99.99%), in 
the sense that massive galaxies of any kind are more metal 
rich than less massive ones. A similar result is found when 
using the mass derived by the fitting to SSP synthetic mod- 
els, instead of the adopted one, based on the Bell & de Jong 
(2001) relation. Although the result is not surprising for the 
E-type galaxies, it is interesting that the L-type galaxies fol- 
low a similar relation, but with an offset. The trend between 
metallicity and mass seems to be universal for any kind of 
galaxy, despite of its morphology. However, for any fixed 
mass the L-type galaxies seem to be more metal rich than 
the E-type ones. This effect is not due to a difference in age 
of the stellar populations, since at equal age we observe a 
considerable range in metallicity. However, our sample is too 
reduced to draw firm conclusions. 



4 DISCUSSION 

The galaxy population of clusters in the local universe are 
domina ted by early-type galaxies ( E 's and bulge-do minated 
spirals) Hubb le fc HumasonI (|l93ll ): lOemleil (|l974l ). with a 
reduced fraction of late-type galaxies (Spirals, Irregulars and 
Mergers). There is a general agreement about the ages of the 
stellar populations in early-type cluster galaxies. It is well 
known that they show a strong color-magnitude relation at 
various redshifts consistents with passive evolution of stel- 
lar populations that formed at z >2-3, with an slope due to 
the mass-metal l icity r e lation (e.g . [Stanford et al.ll 19971 . 1 19981 : 
iKodama et al.1 Il998l ; iPoggiantil |2003| ). The spectrocopic 
features and the scaling relations, like the fundamental- 
plane, the mass-to-light ratio, the Mg&-a and the Faber- 
Jackson relations, explored at different redshift are also 
consistent with passive evolu t ion of an old stellar popula- 
tion (e.g. iBendcr et al.' '1996*; van Dokkum fc Franx ! [l996l: 
Kelson et al. 2001; Zicgler ct al. 2001; van Dokku m et al.l 
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However, this apparently simple scenario for the evo- 
lution of early- type galaxie s in clusters needs to be con- 
trasted with another results. iFasano et ahl l|2000l ) show that 
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there is a substantial morphological evolution in clusters: as 
the redshift decreases the SO population tends to grow at 
the expense of spiral population, whereas the frequency of 
spheroidals remains almost constant. Therefore the popula- 
tion of early-type galaxies observed in distant clusters does 
not nec essarily compris e all the early-type galaxies existing 
at «=0 l|Poggiantill2003l ) . On the other hand, the blue galax- 
ies observed in distant clusters, responsible of the Butcher- 
Oemler effect, must largely become red and fainter. There is 
strong evidence for a quenching of star formation in galax- 
ies as a consequence o f the i nfall in the cluster dense regions 
ijPoggianti et al.ll200ll : IPogg ianti 2004 ), althoug h a previous 
starburst may be also present ( Poggi antil |2004| ) . For galax- 
ies with a star formation history typical of L-type, once 
truncated, it is expected that they become red der and fade 
by 0.5-1.5 mag in ~2-5 Gyr l|Poggiantil [2003). Finally, al- 
though the spectrophotometric and the morphological evo- 
lution seems to be largely decoupled, it is expected that after 
longer timescales there will a morphological t ransformation 
from pure L-type to bulge-dominated spirals l|Fasano et al.l 
l200d ). All these results suggest that the early-type galaxies 
in clusters are composed of a mix of different stellar pop- 
ulations, with very different origins, contrary to the simple 
scenario drawn above. 

How do our results compare with previous knowledge of 
the evolution of galaxies in clusters? The morphological dis- 
tribution of the cluster members derived from our analysis 
is: 11 (7) E-type, 19 (8) lenticulars, and 29 (13) L-type, for 
our complete sample (our IFU sample). In percentages it is: 
19% (25%) E-type, 32% (29%) lenticulars and 49% (46%) L- 
type. These fractions differ strongly with the estimated ones 
for clusters at si milar redshifts: 30% E-type, 50% lenticulars 
and 20% L-type l|Fasano et al.ll2000|). show i ng a m uch larger 
fraction of late- type galaxies. Fasano et al.l l|2000 l) estimated 
these numbers based on images with limiting magnitudes of 
the order of r ~20.2-20.5 mag, i.e. ^^1.5-2 mag less deep 
than our current survey (taking into account the range of 
possible r — z colors for any ki nd of galaxy at this r edshift, 
from ~0 mag to ~0.6 mag, iFukugita et al.lll995l ). If we 
consider only the objects brighter than z <20 mag, the new 
morphological segregation is: 9 (7) E-type, 16 (7) lenticulars 
and 16 (4) L- type. These percen tages are within the values 
reported by (jFasano et al.ll2000l ): 22% (39%) E-type, 39% 
(39%) lenticulars and 39% (22%) L-type. Therefore, there 
is an increase of the fraction of late-type galaxies at fainter 
magnitude regimes. 

Looking at Figures [6] and [12] it is possible to identify 
two kind of L-type galaxies in the faint magnitude regime: 
(i) a group of blue and faint galaxies, with young stellar pop- 
ulations (<8 Gyr), and (ii) another group of red and faint 
galaxies, with older stellar populations (>8 Gyr). The first 
ones are the kind of galaxies responsible for the Butcher- 
Oemler effect. Galaxies recently captured by the cluster, 
with relatively young stellar populations and traces of recent 
star formation processes. Due to their colors these galaxies 
would be detected in surveys in the optical/blue bands. The 
second group are the kind of galaxies predicted in the two- 
steps evolution hypothesis for galaxies captured by clusters, 
where the star formation was quenched first , and then there 
was a slower morphol ogical transformation (|Poggianti et al.l 
I2OOII : |Poggianti|[2003l ). As a first step the late- type galaxies 
fall into the cluster, and due to gas stripping the star for- 



mation process is abruptly interupted (maybe after a short 
starburst process). Then, these galaxies become fainter and 
redder, while still showing a spiral shape. 

Clusters with the highest fraction of blue galaxies shows 
signs of recent mergers, and there are larger fractions of star- 
form ing galaxies in substructured than in relaxed clusters 
fe.g.lSmail et al.lll998l : iMetevier et~aLll2000l : iPimbblet et al.l 
|2002| ) . Theoretical simulations predict that the cluster merg- 
ers induces a starbu rst processes, due to t he strong tidal 
forces involved (e.g IMoss fc Whittle! I2OOOI), a similar ef- 
fect of galaxy-galaxy interactions (e.g. iMihos fc HernquistI 
1996). Abell 2218 sh ows two clear peaks i n the X-ray and 
galaxies distribution (iMcHardv et al.lll990l ) , indicating that 
it most probably has suffered a recent (or ongoing) merg- 
ing. This may account for the moderate Butcher-Oemler 
fraction of blue galaxies quoted in classical studies (~10% 
i Butcher fc Oeml er 1984). 

Galaxies captured by the cluster may suffer similar tidal 
forces, which may induce an enhancement of the star forma- 
tion process, pre vious to the quenching due to gas stripping 
l|Poggiantil |2003| ) . The location of the k-l-a (post-starburst) 
galaxies, and the strength of their Bal mer absorption fea- 
tures strongly supports this hyphotesis l|Poggiantill2004l ) . To 
our knowledge, Abell 2218 is the only cluster at cosmologi- 
cal d istances where HI is observed in emission (|Zwaan et al.l 
l200ll ). at the location of a spiral galaxy in the outer regions 
of the cluster. This galaxy is most probably an infalling spi- 
ral, consuming its gas reservoir with ongoing star forma- 
tion (estimated SFR of <1.4 M0/yr). More spectacular is 
galaxy number 38 from our IFU sample, a blue disrupted 
galaxy at the nominal redshift of the cluster with a very 
little projected distance from the center of the cluster. It is 
possible that this galaxy has been recently captured by the 
cluster, and that it is moving in an orbit perpendicular to 
the line-of-sight. If not, there should be a conspirancy be- 
tween its cosmological redshift and line-of-sight velocity to 
compensate each other so that it appears to be at the nomi- 
nal redshift of the cluster, which is an unlikely explanation. 
Its emission lines indicate a moderate amount of ongoing 
star form ation, with a SFR'~-^6 MQ/yr, derived f rom the Ha 
emission (|Kennicuttl 19831 : iKennicutt et al.lll993 ). No dust is 
derived from the Ha/H/3 ~ ratio, ~2.9, nearly the nominal 
value for case B recombination. Its dominant stellar pop- 
ulation has an age of ~250Myr, which could be a feasible 
time-scale for the infall process. Its morphology is similar to 
that one of a merger in its last stage, however, the strong 
tidal forces that may suffer during the infalling process may 
have produce the same effect. If the proposed scenario for the 
evolution of galaxies in cluster is true, this galaxy would be a 
good candidate for a progenitor of the observed k-l-a galaxies 
at low-redshift (|Poggiantill2004h . In a recent s t udy o f Abell 
2125, a dense cluster at z ~0.247, lOwen et af] l|2006l ) found 
a galaxy with strong emission lines (maybe even with an 
AGN), blue colors, and a disturbed morphology. Although 
that galaxy is brighter than galaxy number 38, both systems 
seem to experience the same physical process. 

Our results indicate that there is a correlation between 
the morphology and the mass, in the sense that less-massive 
galaxies are mostly L-type and more-massive ones are mostly 
E-type (Fig. I13p . This result is natural in the proposed sce- 
nario, where the morphological transformation between L- 
type to bulge-dominated spirals, and later to E-type galax- 
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ies is a consequence of consecutive merging processes. We 
have also found a clear correlation between the mass and 
the metallicity, in the sense that more massive galaxies are 
more metal rich (Fig. I14p . This correlation is an extension 
of the well known Mg-cr relation for spheroidal galaxies, 
which seems to be independent of the galaxy morphology. 
Although we have a very reduced sample, it seems that L- 
type galaxies are more metal-rich than E-type ones at any 
fixed mass. Figure [14] also indicates that the more metal rich 
L-type galaxies are never more metal than the E-type ones, 
being only less massive. To our knowledge this effect has not 
been previously reported in the literature, and it may well 
be an effect of the very reduce sample. If not, it may im- 
ply that the galaxies reach their maximum metallicity when 
they are still spirals, maybe due to the enhancement of the 
star formation rate when infalling into the cluster, and that 
they later evolve to more massive bulge-dominated spirals 
and E-type galaxies keeping this metallicity. 



5 CONCLUSIONS 

We obtained deep IFS data (/h™ ~21.4 mag) covering 
a field-of-view of '^74"x64" centred on the core of the 
galaxy cluster Abell 2218, sampling the wavelength range 
4650-8000A with a resolution of ~10A (FWHM). Combin- 
ing these data with a deep HST/ACS F850LP-band image 
{ziim ~28 mag), and ground-based and HST/WFPC photo- 
metric data, we have obtained the following results; 

• Morphological parameters of all the galaxies in the field- 
of-view of the HST/ACS image down to z <22.5 mag were 
obtained, based on an automatic procedure using galf it. 

• The integrated spectrum of each of the 43 objects de- 
tected in the IFU datacube was obtained, using a modeling 
technique for IFS (galfitSd), which ensures a proper de- 
blending of multiple components in crowded-fields. 

• Of the detected objects, 28 were classified as cluster 
members based on their redshifts. They comprise the first 
volume and fiux limited sample of confirmed cluster mem- 
bers obtained for the core of this cluster. We need to remark 
that the sample was obtained without any pre-selection 
based on morphology, colors or magnitudes. 

• By combining the different previous redshift and photo- 
metric surveys with our new morphological information an 
extended sample of 59 cluster members was created, 31 of 
them sampling the outer regions of the cluster (out of the 
core). 

• The color-magnitude diagram of the full sample of con- 
firmed cluster members, as a function of the morphology, 
was analyzed. We have found that the fraction of early-type 
galaxies is lower than previously reported, most of them 
being red and luminous galaxies. On the other hand, the 
L-type galaxies spread over a wider range of colors and lu- 
minosities with respect to the E-type ones and previously 
published catalogs. 

• There is no significant difference between the Color- 
Magnitude distributions of the galaxies in the core of the 
cluster and in the outer regions if one distinguishes between 
different morphologies. If any, we found a few more L-type 
faint and red galaxies in the core. However, this may be due 
to the biases of the pre-selection of candidates in the samples 
obtained from the literature (in the outer regions). 



• A strongly distorted blue galaxy is found in projection 
at the core of the galaxy cluster. Its redshift corresponds to 
the nominal redshift of the cluster, and therefore, it is more 
likely that it is physically located in the core. A visual in- 
spection of the object indicates that it looks like a merger in 
an advance stage of evolution, dominated by a young stellar 
population. However, if we are observing it in the moment 
of being captured by the cluster, the strong tidal forces that 
it is suffering may produce the observed disruption on its 
morphology. 

• The analysis of the age and metallicity of a reduced sub- 
sample of galaxies (24 galaxies, 12 core members) based on 
line indices indicates that L-type galaxies on average spread 
towards younger ages than E-type ones. On the other hand, 
they seem to cover the same range of metallicities at any 
fixed age. 

• The age and metallicity of each galaxies in the core 
of the cluster was derived by a fitting procedure of its in- 
tegrated spectrum to single stellar population models. The 
age- metallicity diagram of these galaxies indicates too that 
L-type galaxies cover a clearly wider range of ages, and a 
slightly wider range of metallicities. The more massive L- 
type galaxies seem to have similar ages and metallicities 
than the E-type ones, while the less massive ones tend to be 
younger and less metal rich. 

• The distribution of ages and Sersic indices as a function 
of the stellar masses of the core members shows that L-type 
galaxies are less massive on average. They spread over a 
wider range of ages, but it is not possible to determine if this 
is due to their morphologies or masses. On the other hand, 
more massive galaxies tend to have older stellar populations 
and they are all early-type galaxies. 

• There is a clear trend between the metallicity and the 
mass, in the sense that more massive galaxies are more metal 
rich. This trend is just a change of coordinates of the well 
known Mgfc-cr relation for the early-type galaxies. However, 
we found that this trend is also present for the L-type galax- 
ies, but with an offset: they are more metal rich than E-type 
galaxies at any fixed mass. 

Our results agree with the proposed t wo-steps scenari o 
for the evolution of galaxies in clusters (e.g. |Poggianti|[2"003l ') . 
in which we have (i) a primordial population of early-type 
galaxies formed at early epochs with stellar populations that 
are mostly old, and (ii) a new population of L-type galaxies 
that are captured by the cluster, infalling, suffering a short 
enhancement of the star formation that is later quenched 
by the interaction with the environment. The galaxies then 
evolve passively, becoming redder and fainter, while keep- 
ing the spiral morphology. The disks dim due to the lack of 
new formed stars, and the galaxies become bulge-dominated 
spirals, after longer time-scales. Later galaxy-galaxy inter- 
actions and dry merging processes build up new massive 
spheroidal galaxies. 
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APPENDIX A: INDIVIDUAL SPECTRUM 

Figure lAll shows the indivual spectrum extracted from 
our IFS data for each galaxies listed in the Table |4] For 
each galaxy the spectrum is plot together with the multi- 
band broad-band photometry derived from the HST/ACS 
F475W, F555W, F625W, F775W and F850LP-band images 
(red circles). The vertical error bars indicate the photomet- 
ric errors, while the horizontal ones indicate the width of 
each band. The wavelength of the strongest night sky emis- 
sion lines and bands are also indicated. At these locations 
the unproper sky subtraction produce characteristic residu- 
als that can be seen in the spectra. In particular the wave- 
length range redder than 7200 A is strongly affected, so we 
did not consider it in any of the performed analysis. The 
grey solid line shows the best fitted SSP model derived with 
the procedure described in Section 13.7.21 
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APPENDIX B: MORPHOLOGICAL ANALYSIS 
AND SPECTRA EXTRACTION: CHROMATIC 
EFFECTS 

The procedure to deblend and extract the spectra of the in- 
dividual objects from our IFS data relies on the assumption 
that it is possible to use the morphological parameters de- 
rived by fitting galaxy profile models on the F850LP-band 
image, using galf it, to derive the photometry at any other 
wavelength range, by forcing galf it to recover the flux with 
the morphological parameters fixed. To test this assumption 
use the HST/ACS multiband image dataset, which cover all 
the wavelength range or our IFS data. 

For each broad-band image of our dataset (F475W, 
F555W, F625W and F775W), we performed the same 
morphological analysis that we already performed for the 
F850LP-band, following the prescriptions described in Sec- 
tion 13.21 In this first iteration we left galf it to fit the 
morphological parameters, as we did in the analysis of the 
F850LP-band. This procedure produces a morphological 
catalogue for each of the considered images, including the 
same parameters listed in Table [S] In a second iteration, the 
analysis was repeated but this time we used as input parame- 
ters the ones derived by the analysis of the F850LP-band im- 
age. In this second step the morphological parameters were 
fixed to those input values, and only the photometry was 
fitted by galf it. This procedure is basically the same that 
we use to extract and deblend the spectra in our IFS data. 

Similar results have been found for any band. Figure [Bll 
summarizes the results from the analysis. Left-panel shows a 
comparison between the different bands magnitudes derived 
by galf it forcing the morphological parameters to be the 
ones derived from the morphological analysis of the F850LP- 
band image, and the magnitudes derived not imposing any 
condition to galf it regarding the morphological parameters 
of the analyzed objects. The magnitudes derived using both 
methods agree within ±0.2 mags, without any clear trend or 
bias. Middle-panel shows the scale-length derived by galf it 
for the F850LP-band image, using in the forced fitting pro- 
cess over the different bands' images descrived above, along 
the scale-length derived directly over this image by the mor- 
phological analysis. Although there are some outlayers, both 
scale-lengths are similar in most of the cases. Right-panels 
shows a similar comparison between the Sersic indices de- 
rived by the two fitting procedures. In this last case we found 
the biggest differences between both parameters, with a con- 
siderable dispersion from the one-to-one relation. It is clear 
that there is a chromatic effect over the morphology, as it 
is seen in the last panel, however, this effect has a reduced 
impact on the derived scale-length, and even less impact on 
the derived photometry. Thus, the basic assumption of our 
procedure to deblend and extract the individual spectra of 
the targets from our IFS data is valid. 

To perform a final test of the validity of our pro- 
cedure and the accuracy of our spectrophotometric data 
we synthetized the fiuxes throughout the B, V, R and I- 
bands for each of our 28 cluster members with IFS data by 
convolving their extracted spectra with the corresponding 
HST/ACS filter band transmission curve, integrating their 
fluxes. These fluxes are compared to the corresponding ones 
derived using the HST/ACS based photometry listed in Ta- 
ble m Figure IB2I shows the histograms of the relative differ- 



ences between both fluxes for each band, AF/F, that basi- 
cally corresponds to the error in the corresponding magni- 
tude (~emag). Due to the high accuracy of the HST/ACS 
photometry (typical Cmag ~0.05 mag for the considered ob- 
jects) most of the differences shown in the figure corresponds 
to spectrophotometric errors in the IFS data. For all the 
bands both fluxes agree within a range of ±0.3 mag, with a 
standard deviation of the difference of ~0.15 mag. 
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Figure Al. Individual extracted spectrum for each of the confirmed cluster members within the field-of-view of our IFS data (black 
solid line) , each one labeled with the identification number listed in Table |4] For each galaxy the spectrum is plot together with the 
multiband broad-band photometry derived from the HST/ACS F475W, F555W, F625W, F775W and F850LP-band images (red circles). 
The vertical error bars indicate the photometric errors, while the horizontal ones indicate the width of each band, he grey solid line 
shows the best fitted SSP model derived with the procedure described in Section l3.7.2l 



The galaxies in the core of Ahell 2218 25 




4000 5000 6000 7000 8000 9000 

Wavelength (A) 



A2218 11 oky Em. Lines Sky Lm bands 



"1 I I 



1 




4000 5000 



6000 7000 8000 

Wavelength (A) 



4000 5000 



6000 7000 8000 

Wavelength (A) 



I 
S 



T ' r 

1 ^T^T 



(oaD nma tnao ivx 



4000 5000 6000 7000 

Wavelength (A) 



8000 9000 




A221B 12 any Em- Lines Sky Em Bands 

^1 I I I I I 




9000 4000 6000 6000 7000 8000 9000 

Wavelength (A) 



A221B 13 Sky ^m. Lines Sky Em. Bands 

r 




lA - 



-r 





ASIA Ifl r^y 




4011C 31HX1 HOOD TWn BODC BDDD 

WbViiLf-D^lL [ill 



.——I . 1—. — r- 

B.: Un>> Shr Itcn tl«ni'l* 



4000 ^[HXl 



DODD TWn 
WEMelciijlli 1!LJ 
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Figure Bl. left-panesl: F475W, F555W, F625W and F775W-band magnitudes derived by galf it forcing the morphological parameters 
to be the ones derived from the morphological analysis of the F850LP-band image, as a function of the same magnitudes derived not 

imposing any condition to galf it regarding the morphological parameters of the analyzed objects, middle-panels: similar comparison for 
the length-scales derived by galf it for the F850LP-band image and for any other band . right-panels: similar comparison for the Sersic 
indices. 
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Figure B2. Relative differences between the fluxes obtained for the B,V ,R and /-bands derived from the IFS extracted spectra and 
those ones derived using the photometry obtained from the HST/ACS images. 



